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Abstract. Much of our knowledge about the solar dynamo is based on sunspot observations. It is thus 
desirable to extend the set of positional and morphological data of sunspots into the past. Gustav Sporer 
observed in Germany from Anklam (1861-1873) and Potsdam (1874—1894). He left detailed prints of 

■ sunspot groups, which we digitized and processed to mitigate artifacts left in the print by the passage of 
" time. After careful geometrical correction, the sunspot data are now available as synoptic charts for almost 

450 solar rotation periods. Individual sunspot positions can thus be precisely determined and spot areas can 

■ be accurately measured using morphological image processing techniques. These methods also allow us to 
1 determine tilt angles of active regions (Joy's law) and to assess the complexity of an active region. 
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Q_[ The first continuous observations of sunspots were made by Samuel Heinrich Schwabe 1825— 

1867. In this time, Richard Carrington started his observations in 1853, followed by Sporer's 
$_i ■ first observations in January 1861 during Carrington rotation period No. 96. In total 445 rotation 
c/3 . periods were recorded. The images were published as lithographs, which reproduce many details 
of the sunspot fine structure, e.g., umbra and penumbra (Sporer 1874). Each page contains five 
rotation periods. The abscissa shows the heliographic longitude L from 360° to 0° and on the 
ordinate is the heliographic latitude b from +40° to —40° for one rotation period. 

The sunspot drawings were collated in hardbound books. Consequently, pages will not be 
perfectly flat when opening the books. Binder clips and metal coins were used to get a flat im- 
age of the drawings. Whole pages were photographed with a digital camera. The drawings were 
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1. Sunspot observations and digitization of sunspot drawings 



■ published nearly 150 years ago, thus, the pages are yellowish, stained, and corrugated. These ar- 



tifacts have to be removed because they would severely affect the analysis of sunspot properties. 
All algorithms were developed in the Interactive Data Language (IDL). 

2. Image processing and sunspot morphology 

After converting the color images to gray-scale images, they are rotated by 90° to have them 
/\ \ in portrait format as in the books. The images are cropped so that only the coordinate grids with 
£2 ■ sunspot data and their annotation remain. To remove a small rotation angle error of the camera, 
the standard deviation of the intensity traces of all rows and all columns has to reach a maximum. 
Then, the angular offset is zero. An order-statistics filter (Gonzalez & Woods 2002) is used to 
compute a background intensity map. The filter ranks the intensities within a 33 x 33 pixel 
neighborhood and replaces the value of the center pixel by the 80 th percentile of the intensity 
distribution. The 16-pixel wide borders of the full background intensity map are replaced by the 
mean background intensity. Division of an image by its background intensity map (i.e., a pseudo 
flat-fielding procedure) results in an artificially bleached and flattened image. 
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Figure 1. Image after remapping, where crosses mark the sunspot positions and circles the mean positions 
of the sunspot groups. The positions were measured over several days, i.e., it is possible to visualize the disk 
passage of sunspot groups. Sometimes there are errors in the data, e.g., the table entry of spot 7 differs from 
its position in the image. The insert shows the results of blob analysis for active region 14. 

To get the exact coordinates of the grid lines, a 17 x 17-pixel kernel is used, where the value of 
two intersections will be set to unity and all other pixels at zero. A good approximation of the grid 
coordinates is the average of the line profiles of rows and columns, which were smoothed with 
a 17-pixel Lee-filter and afterwards normalized. Good estimates of the x- and y-grid coordinates 
yields a parabolic fit to local maxima of row and column profile. To improve the grid coordinates, 
from the binary mask, a region of 20 pixels will be extracted and enlarged by factor of five. 
Afterwards, the region will be smoothed with a 5-pixel boxcar filter. Nevertheless, sunspots near 
to intersecting grid lines can lead to erroneous results. As a precaution, all x- and ^-coordinates 
are checked against the positions of their immediate neighbors and replaced, if a significant offset 
is detected. The heliographic longitude and latitude are given in increments of 10°, so that we 
have a grid of (37 x 9) x- and y-coordinates for each rotation period. We also included a 5°-wide 
boundary to minimize sampling errors at the periphery. Subsequently, Delaunay triangulations 
were used to interpolate the sunspot data to the new grid (the result is shown in Fig. 1). 

After image processing, standardized sunspot maps are at hand for further analysis. The sunspot 
drawings are made of dark-inked areas for the umbrae and many single points to indicate the 
penumbral fine structure. A more realistic 'intensity distribution' can be achieved using some 
Gaussian smoothing. Even though blurred, the images are now suitable for image processing 
techniques such as intensity thresholding. Morphological opening and closing operations can 
then be applied to the binary masks, in which individual sunspots are labeled. Thus, 'holes' in the 
sunspot area can be filled and spot contours become less rough. We used standard tools for 'blob 
analysis' (Fanning 201 1) to derive parameters describing the morphology of sunspots and to de- 
termine their positions. In general, this approach leads to very good results but in some cases, the 
spoke-like structure used to indicate penumbral fine structure results in artifacts. Another prob- 
lem are the numbers, which label the sunspot groups. However, since they consist of thin lines 
as compared to the more compact sunspots, morphological image processing is again capable of 
indentifying and removing the labels. In summary, Sporer's sunspot drawing are of high quality 
facilitating the application of modern feature recognition and image processing methods. 
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